AD020945 


APPLICATIONS  OF  ANALYTIC  CONTINUATION  TO  THE  SOLUTION  OF 
BOUNDARY  VALUE  PROBLEMS 


BY 

P.  R.  GARABEDIAN 


TECHNICAL  REPORT  NO . 16 
NOVEMBER  20,  1953 

/ 

PREPARED  UNDER  CONTRACT  Nonr-225(ll) 

(NR- 041- 086) 

FOR 

OFFICE  OF  NAVAL  RESEARCH 


APPLIED  MATHEMATICS  AND  STATISTICS  LABORATORY 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA 


APPLICATIONS  OF  ANALYTIC  CONTINUATION  TO  THE  SOLUTION  OF 


BOUNDARY  VALUE  PROBLEMS 
By 

P.  R.  Garabedian 

A better  analysis  of  special  solutions  of  particular  boundary  value 
problems  is  of  considerable  importance  in  the  theory  of  partial  differential 
equations,  especially  for  equations  of  fourth  order,  because  very  few 
explicit  solutions  are  known  in  this  case.  The  present  work  was  motivated 
by  interest  in  the  non-linear  Navier-Stokes  equation 

zzzz  z zzz  ^z  zzz 

governing  the  steady  two-dimensional  flow  of  an  incompressible  viscous 
fluid,  but  the  results  have  application  in  other  fields,  such  as  elasticity. 
We  develop  methods  based  on  analytic  continuation  which  yield  for  certain 
special  plane  domains  the  solution  of  the  basic  boundary  value  problems 
for  the  equations  A<p  ■ qp  and  A A ^ ■ 0.  The  methods,  which  consist  either 

•R> 

in  the  introduction  of  independent  complex  variables  z-x+iy  and  z - x-  iy 
or  in  a suitable  exploitation  of  the  analytic  function  G(z)  such  that  the 
equation  z«G(z)  represents  the  boundary  of  the  domain,  are  of  themselves 
interesting,  since  they  provide  formulas  for  the  reflection  of  solutions 
of  these  partial  differential  equations  across  analytic  boundaries.  We  are 
thus  able,  for  example,  to  study  the  solution  of  the  Oseen  equation 
near  the  end  of  a slit. 

Let  P(z,z)  be  a fixed  real  analytic  function  of  x and  y,  and  let  cp(z,z 
be  a real  solution  of  the  linear  elliptic  partial  differential  equation 


-W — 
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(l)  Cp^_(z,z)  - P(z ,z) cp(z,z) 

in  a plane  domain  D,  where 

_g2£  _a_  . l(.a_  _ t .*_)  , m-2L)  . 
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We  denote  by  A(z,z  jt,t  ) the  Riemann  function  of  the  corresponding  hyperbolic 


equation 

(2) 


d2q> 


- P(z,z*)(p 


dz  dz 

# 

:Ln  two  independent  complex  variables  z and  z . In  other  words,  A is,  as 
a function  of  z and  z* , a solution  of  (2)  which  satisfies  the  boundary 
conditions  A(z,t#;t,t  ) - A(t,z*jt,t*)  ■ 1.  If  we  continue  the  solution 
q?(z,?)  of  (l)  analytically  into  the  domain  of  complex  values  of  x and  y, 
it  becomes  a solution  Cp(z,z  ) of  (2).  It  is  clear  by  Stokes'  theorem 
that  because  <p  and  A solve  (2),  the  lino  integral 


(3) 


/{* 


dz 


dz 


dz*J- 


is  independent  of  the  path  of  integration  between  given  limits.  We  make  use 
of  this  result  to  derive  a dual  pair  of  useful  formulas  expressing  <p  in 
terms  of  A. 


The  required  formulas  are 

t 


- . ,*x 


(4)  <p(t,tW)  - cp(t,t)  -j\p(z,z)  1 ds+A(z,z;t,t  ) dz^-  , 


_# 


where  the  path  of  integration  runs  from  t to  t in  the  real  plane  domain  D, 
and  [1] 

t * 

(5)  (p(t;t)~  ^>\t  jt  )A(t  ,t  jt,t )*J  A(z,t  jt,t)  ^ dz 

t*  € 

+ f A(t*,z*jt,t)  dz*  . 


t 


a z 
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The  first  formula  is  derived  by  evaluating  the  integral  (3)  over  a closed 

* _ -* 

circuit  consisting  of  a curve  in  the  real  plane  z - z from  t to  t,  plus  a 

* _ * * 

curve  in  the  characteristic  plane  z « t from  z *t  to  z - t , plus  a curve 

* * -* 

in  the  characteristic  plane  z • t from  z - t to  z * t , for  the  integrals  in 
the  characteristic  planes  simplify  by  virtue  of  the  boundary  conditions 
satisfied  by  A.  In  a similar  way,  formula  (5)  follows  by  evaluation  of 

(3)  over  a closed  path  made  up  of  the  four  arcs  which  .loin,  respectively, 

* - * * -* 

the  points  z «t  and  z - t in  the  plane  z*t,  the  points  z • t and  z * t 

* * * * * _ -* 
in  the  plane  z ■ t , the  points  z - 1 and  z ■ t in  the  plane  z * t , and 

the  points  z-t  and  z-t  in  the  plane  z*«t.  Finally,  since  (f>( z,i) 

and  A(z,z;t,t)  are  real,  (5)  reduces  to 

t * 

(6)  Cp(t,t)  - <p( t ,t  )A(t  ,t  ;t,t)  ♦ 2 Re  -f  J A(z,t  ;t,t)  dz'j'. 

t* 

As  a first  application  of  the  formulas  (4)  and  (6)  we  consider  the 
equation 

(7)  Acp-cp 

in  the  wedge  - ctf  < arg  z < oC  . In  this  case  P - l/4  and  the  Rieraann 
function  A has  the  form 

A(z,z*;t,t  ) - J[i^(s-t)(z*-t*)'  ] , 

<|fr 

where  J is  the  Bessel  function  of  order  zero.  We  set  t -0  and  we  take  as 
path  of  integration  the  ray  from  the  origin  to  z,  whence  (4)  yields  for 
the  solution  ( p of  (7) 

t 

(8)  <^(t,0)  - <^(t,t)  -j  i 9‘^z  ]dz+  J[  t-z)z  * ] dz\  , 
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The  remarkable  fact  about  this  special  case  is  that 


t 

(9)  He{<p(t,0)^  - <f(t,t)-  + | d'z'j 

^ 0 

since  (t-z)z  >0  along  our  path  of  integration.  Thus,  along  any  ray  through 
the  origin,  Re-|9(t,o)^  is  completely  determined  by  ^>(t,t)  alone.  Further- 
more, for  real  arguments  the  Bessel  function  J is  uniformly  bounded, 
whereas  Cp(t,t)  and  its  derivatives  are  required  to  vanish  exponentially  at 
infinity  for  solutions  of  boundary  value  problems  in  a wedge.  Therefore,  by 
(8),  <^(t,0)  will  be  bounded  for  such  solutions.  It  follows  from  this 
analysis  that  the  first  boundary  value  problem  for  (7)  in  a wedge 
- ct  -c  arg  z transforms  by  (9)  into  a Dirichlet  problem  for  the  harmonic 
function  Re|cp(t,0)^j'  in  the  same  wedge.  Thus  <^>(t,0)  can  be  found  by 
standard  procedures  and  the  explicit  solution  <^t,t)  of  (7)  assuming  the 

given  boundary  values  can  be  obtained  as  an  integral  from  the  LOirwlila 

t 

(10)  cp(t,t)  - cp(0,0) J(i  1 1 1 ) + 2 Rs  JCl  a]  (t-z)t*  ] dz  , 


based  on  i6) . 

The  general  formula  (6)  for  solutions  Cp  of  (l)  in  a domain  D gives 
a convenient  procedure  [6]  for  continuing  Cj>  analytically  across  an 
analytic  arc  of  the  boundary  of  D along  which  Cj?  satisfies,  for  example, 
the  boundary  condition  cp-0.  We  write  the  equation  of  the  analytic  arc  in 
the  form  z»G(z),  where  G(z)  is  an  analytic  function.  Such  a representation 
can  always  be  obtained  by  solving  for  z the  more  usual  equation  of  the  arc. 
We  suppose,  without  loss  of  generality,  that  the  origin  is  the  point  on 
the  boundary  of  D near  which  we  wish  to  make  the  reflection.  In  D,  near  the 
origin,  consider  the  analytic  function 
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(11) 


V 

0(t ) - J A(z,0;tfG(t))  dz 

o 

From  (6)  and  the  boundary  condition  Cp  * C,  we  obtain  along  the  curve  z*G(z) 
the  identity 

. t 

0 - qpU,t)«  <p(0,0)A(0,0jt,t)  + 2 Re  -j  A(z,Ojt,t)  ^ — dzj 

- 2 Reju(t)^  -D(t)  + U(t) 

Hence  the  analytic  function  U(t)  can  be  continued  analytically  across  the 
curve  z~-  G(z)  by  the  rule 

U(t)  - -U(G(t) ) , 

since  when  t lies  on  one  side  of  the  curve  and  sufficiently  near  the  origin, 
the  point  G(t)  lies  on  the  opposite  side.  We  now  derive  from  (ll)  the 
integral  equation 

t 

(12)  0(G(t))*jAU,0jt,G<t))  9jf^0)  dz  - 0 

o 

for  the  determination  of  the  analytic  function  <p(t ,0)  across  the  arc 

z«G(z).  This  equation  is  linear  and  of  the  Volterra  type,  so  there  is 

no  difficulty  is  solving  it  by  successive  approximations.  Thus  the  analytic 

continuation  of  Cf>( t,0)  is  performed,  and  tc  obtain  the  continuatici/of 

the  original  solution  Cjf>( t,t)  of  (l)  we  have  only  to  apply  formula  (6) 

# 

again  with  t - 0. 

The  procedure  just  outlined  can  readily  be  extended  to  the  case  of 
fourth  order  elliptic  equations . We  illustrate  this  generalization  with 
an  application  to  the  important  special  case  of  the  Oseen  equation 

(13) 

This  equation  has  a general  solution 
^ * ex  cp+  h 


y 
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where  qp  is  a solution  of  (7)  and  h is  a harmonic  function.  By  (10),  the 
general  solution  can  be  written  near  the  origin  in  the  form 


t 

(14)  " Re  J[i  ^ (t-z)t1  ]g'(z)dz  ♦ Ce^  ^’^2J(i|t  |)  + f(t)^ 


o 

where  C is  an  arbitrary  constant  and  f(t)  and  g(t)  are  arbitrary  analytic 
functions.  We  suppose  that  the  solution  4 (13 ) is  defined  in  a 

region  D whose  boundary  contains  the  analytic  arc  z-G(z)  passing  through 
the  origin,  and  along  this  arc  we  assume  that  the  boundary  conditions 


(15)  4*  - - 0 

do 

are  fulfilled,  where  n denotes  the  inner  normal  of  the  arc.  The  problem 
is  to  continue  4 analytically  across  the  arc  z-G(z),  and  it  is  clearly 
sufficient  by  (14)  to  perform  such  a continuation  on  the  analytic  functions 
f (t)  and  g(t) . 

We  denote  by  s the  arc  length  along  the  curve  z-G(z),  and  we  verify 

that 

(16)  ’ 1 - iG'(z)"l/2 

ite placing  ( 15 ) by 

e”x  4-  die~xy)/-dn  - 0 
and  substituting  (14)  into  the  result,  we  obtain 


(17)  Re  £v(t) ^ r*  Re £w(t)^  - 0 


along  the  curve  z*G(z),  where  V and  W are  the  analytic  functions  defined  by 


(18)  V(t) 


J J[i  /')(t-z)G(t)l]g','(z)dz  + CJ(i^ tG(t)#  ) + f(t)e  ^^(t))^ 


l 


9 
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(19)  W/(t)G/(t)"1^2  - iG/(t)"1^2|g/(t)  + iCJ'(i  (tG(t)1) 

L 2/\[tG(tT 

+ i fj/(i^(t-z)G7t)1  ) g/(z)dz 

</  2 (t-z)G(t) 

o 1 

in  D.  Therefore  these  functions  can  be  reflected  across  z-G(z)  by  the  rule 
V(t)  - -V(G(t))  , W(t)  - -I(G(t)) 

It  follows  then  from  (18)  and  (19)  that  f(t)  and  g(t)  can  be  continued 
analytically  across  z-G(z)  as  solutions  of  the  system  of  Vo  1 terra  integral 
equations 

(20)  iG/(t)»/(G(t))-g/(t)  + e“^t+G(t))/2f/(t)  + f(t)e'(t+G^t)^2 

♦if  J'U^t-zWt)')  /(z)dz 

i 2 ^ (t-z)G(t)’ 

♦ iCJy(iV  tG(t) 1 ) , 

’ 2*ftG(tT 

(21)  ^CV/(G(t))  + fl'(G(t))]  - e“^t+G^t^/2f/(t)-f(t)e‘^t+G^t^/2  — 

ax  2 

- 1 a f r(^(fc»j?(ii>,><E,/(.)dl.  10  a LiiMjhR . 

dt  J ^G(t)  dt  ^fG(t) 

Once  .the  solution  of  (20)  and  (2l)  has  been  obtained  by  successive  approxi- 
mations, the  continuation  of  the  stream  function  4^  can  be  found  from  the 
express$>n  (14) . | jit 

We  can  apply  the  rule  just  described  for  reflecting  4^  to  the  case 
where  4^  is  defined  in  a neighborhood  of  the  origin  slit  along  the  positive 
real  axis,  on  which  the  boundary  conditions  (15)  are  imposed.  Here  G(z)»z 
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end  it  is  found  that  V and  W are  regular,  single-valued  functions  of /ft  near 
the  origin,  whence  by  (20),  (2l)  the  same  is  true  for  f and  g.  Therefore 
is  a regular  function  of  Nt  and  and  since  ^ and  oust  vanish 

* y 

at  the  origin,  ws  conclude  that  the  skin  friction  AH'  behaves  like  l/'fx 
along  the  positive  real  axis  C 5 1 . 

Our  mothod  of  reflection  is  valid  for  the  equation  AA'H'**  A^H'  of 
the  vibrating  clamped  plate,  also,  and,  indeed,  many  further,  more  involved 
applications  could  be  carried  out. 

We  turn  our  attention  to  the  simplest  fourth  order  equation*  namely, 

(22)  AAH^-  0 , 

and  we  show  that  the  above  technique  simplifies  to  such  an  extent  that  we  are 
led  to  ths  explicit  solution  of  specific  boundary  value  problems . The 
general  solution  of  the  biharmonio  equation  (22)  has  the  form 

(23)  H-’  - Re  ^zf(z)  + g(z)|  , 

where  f and  g are  arbitrary  analytic  functions  of  z.  If  we  substitute  this 
formula  into  the  boundary  conditions  (15)  along  the  arc  z-G(z)  and  take 
(16)  into  account,  we  find 

(24)  Re-^G(z)f(z) + g(z)^r  *0  > 

(25)  Re^iG(z)f/(z)G/(z)"1//2+ ig^(z)G/(z)"1//2-if(z)G/(z)1/,2j' - 0 

We  define  the  two  analytic  functions 

$(z)  - G(z)f (z)  + g(z)  , 

’J'(z)  - iJ,^G(z)f  '(z)  + |'(j5)-f(z)G'(zfy  dz  i, 
and  we  note  that  (24)  and  (25)  imply 


(26) 


Re  ^(z)y  - Re  |^(z)^-  = 0 
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ou  the  curve  s-Gfjs).  Bui  also 

f(z)  - [<$'(*) + i$/(*)]/(2G'(»)) 

g(z)  - $(z)-C<^/(z)  + ilJ,'(z)lG(z)/(2G/(z))  , 

whence  the  biharmonic  function  4 satisfying  the  boundary  conditions  (15) 
can  be  expressed  in  the  fora 

(27)  ^ - Re  {$(«)  + [$'(»)  4 1 $*(»)]% 

in  terms  of  the  pair  of  analytic  functions  $ and  '4  satisfying  the 
equivalent  boundary  conditions  (26). 

As  in  our  earlier  studies,  the  formula  (27)  gives  a rule  for  continuing 
4*  analytically  across  an  arbitrary  analytic  arc  z:-G(z).  Indeed,  we  have 
only  to  reflect  and  Q by  the  usual  Schwarz  principle,  using  (26). 

However,  in  the  present  situation,  the  reflection  rule  is  so  simple  and 
elegant  that  for  certain  domains  D it  can  be  exploited  in  order  to  solve 
explicitly  ths  first  boundary  value  problem  for  (22).  This  can  be  done 
when  D is  bounded  by  a simple  closed  analytic  curve  z-G(z)  such  that  G(z) 
is  single- valued  in  D and  regular  there  except  for  a finite  number  of  poles. 
The  functions  $ and  4^  defined  by  (27)  re  regular  in  D except  for 
singularities  at  the  poles  of  G,  and  j the  boundary  curve  we  find 

(28)  4-Re££} 

(29)  /■ !r as  - MX)-  • 

Hence  the  problem  of  determining  a solution  of  (22)  in  P with  prescribed 
values  for  4*  and  9 4/3  n on  the  curve  bounding  D reduces  to  the  problem 
of  finding  analytic  functions  and  4 with  given  real  parts  on  the 
boundary  of  D end  with  suitable  singularities  at  the  poles  of  G.  The 
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solution  therefore  involves  only  the  determination  of  a finite  number  of 

parameters  associated  with  the  singularities  in  such  a way  that  <$/  ♦ i'J' 

vanishes  at  the  zeros  of  G'.  In  particular,  if  G has  only  simple  poles, 

then  has  only  simple  poles  and  $/+  i<^'  has  only  simple  poles. 

As  an  example  of  this  theory,  we  develop  a representation  for  the 

biharmonic  Green's  function  P(z,^)  of  the  exterior  of  the  ellipse 
2-2  2-2 

xch  od  + ysh  oC  -1.  The  result  should  be  compared  with  the  earlier 
solution  given  by  a related  method  [7,8].  We  write  the  equation  of  the 
ellipse  in  the  complex  form 

(z  + z)^ch”^oC  - (z - z)  sh  = 4 


and  solve  for  z as  a function  of  z to  obtain 
(30)  G(z)  ■ z ch  2eC  - sh  2oc 


Denote  by  p(z,^)  the  analytic  function  of  z whose  real  part  is  the  Green's 
function  for  Laplace's  equation  in  the  exterior  of  the  ellipse,  with  source 
point  at  z*  A fundamental  solution  of  (22)  is 
Re  |z-  ^|2p(z,  , 

and  we  must  add  to  this  an  expression  of  the  type  (27),  adjusted  so  that  the 
sum  satisfies  the  homogeneous  boundary  conditions  (15)  on  the  ellipse.  This 
gives  on  the  ellipse 

Re{£}  - 0 

Re'0^]r  ■ **{^1  1 ^ >dz^  , 

whence  wo  have 


(31)  A(z2-G(z)2)  , 

'£■  i|(z-^)(G(z)-^)p'(z,^)dz  + B(z-G(z))+iC(z+G(z))+Dp(z,oo)  , 


(32) 


11  - 


where  the  real  parameters  A,  B,  C,  D are  to  be  found  from  the  condition 
that  ^5/  + iJ'  must  vanish  at  the  zeros  ±ch2p(.  of  G'U).  We  thus  obtain 
for  A,  B,  C,  D the  equations 

+ 2A  ch  2oC  i i£p'(ch  2o£,  oo)  + i(B+  iC) 


- (ch  2ot  X Xf  )(l+  Y.  )p '(+  eh  2 od . k) 

and  therefore 


(33)  T(z,^)  > Re{-|z-^|2p(z,$)  + $ + g7^j  C^'  + ii"]} 

with  G,  $,  and  £ given  by  (30),  (3l),  and  (32). 

The  above  result,  which  is  of  some  interest  for  the  discussion  of 
slow  viscous  flow  around  an  ellipse,  simplifies  considerably  in  the  limiting 
case  when  oC  - 0 and  the  ellipse  degenerates  into  a slit  from  -1  to  +1.  Here 
we  find,  in  fact,  that  the  biharmonic  Green's  function  is  given  by 

(34)  fU,^)  - Re  -£-|  z-Xf  |2p(z,^)  + iy(z-  ^)(z-  ^)p'(z,  l^)^  , 


since  G(z)«  z.  This  same  formula  yields  the  biharmonic  Green's  function 
for  the  exterior  of  any  finite  number  of  slits  along  the  real  axis,  provided 
that  we  interpret  p to  be  that  analytic  function  whose  real  part  is  the 
harmonic  Green's  function  of  the  multiply- connected  slit  domain. 

With  our  present  method  we  can  treat  the  first  boundary  value  problem 
for  (22)  in  domains  bounded  by  slits  along  an  analytic  arc  z=G(z). 

However,  for  the  applications  it  is  also  desirable  to  study  the  behavior  of 
biharmonic  functions  near  a point  of  the  boundary  where  two  analytic  arcs 
intersect  at  an  arbitrary  angle.  In  order  to  obtain  n first  analysis  of 
problems  of  this  type,  we  determine  here  the  biharmonic  Green's  function 
Hz,^)  for  an  arbitrary  crescent  domain  bounded  by  two  circular  arcs. 


Without  loss  of  generality,  we  can  assume  that  the  two  circular  arcs 
bounding  the  crescent  intersect  at  the  points  1 and  -1.  The  conformal 
transformation  z * thw  therefore  maps  the  crescent  in  the  z- plane  onto  an 
infinite  strip  a < v b in  the  w-plane,  w-  u + iv.  The  general  solution 
(23)  of  (22)  can  be  replaced  in  the  w-plane  by  an  expression 

(35)  F(w,w)  - (ch  w ch  i)^  ■ Re -^f  (w)sh  w + g (w)chw*y  , 


which  is  easily  seen  to  be  the  general  solution  of  the  transformed 
differential  equation 

(36)  (-^  - l)(“o  - l)F  - 0 

The  Green's  function  P(z,  tf)  of  (22)  for  a crescent  will  be  related  to  the 
Green's  function  G(w,(J")  of  (36)  for  the  corresponding  strip  by  the  identity 


(37) 


G(w,  (T)  - I ch  w |2|ch  (T|2  P(z,^) 


where  * th  <T  , (T  ™ s ♦ it.  We  proceed  to  determine  G by  the  method  of 
Fourier  transforms,  omitting  details  of  a familiar  nature  [2,4,91. 

We  write  G(w,(J”)  as  the  Fourier  transform  with  respect  to  0 of  a 
quantity  #(v,t;0),  which  turns  out  to  be  the  Green's  function  with 
discontinuity  at  v*>t  of  the  ordinary  fourth  order  differential  equation 

(38) 

in  the  interval  a -c.v-^b,  since  equation  (36)  has  constant  coefficients. 
Thus  we  have 

CD 

(39)  G(*,<3")  m yf  J ^(v,t;0)cos0  (u-s)dQ  , 

o 

where  the  symmetric  Green's  function  ^(v,t;0)-  ^f(t ,v;0)  is  given  for 


+ (8  - 202)  ♦ (16  + 802  + 04)  - 0 
dv^  dv^ 


a^v^t^b  by  the  lengthy  relation 
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(40)  4©(©2  + 4)(4sh 2 £9-  e2sin22^)  ^(v,t;6) 

- (e2+4)(ech£esin2i-2  ah^e  cos  2i,)[sh©(v-a)sin  2(v-a)sh©(t-b)sin  2(t-b)] 

+ (©2+4)sbiesin  2i’^sh©(v-a)sin2(v-a)[eche(t-b)sin2(t-b)-2  sh©(t-b)cos  2(t-b)] 

- sh©(t-b)sin  2(t-b)(©ch  ©(v-a)sin2(v-a)-2  sh©(v-a)cos  2(v-a) 

- (©chdesin  2j£  +2  sh^9  cos  2 0)[©  ch©(v-a)sin  2(v-a)-2sh  ©(v-a)cos  2(v-a)  1* 

•C©ch  0(t-b)sln2(t-b)-2  sh©(t-b)cos  2(t~b)  ] , 

in  which  we  have  replaced  b~a  by  . 

The  explicit  representation  of  the  biharmonic  Green's  function  for  a 
crescent  given  by  formulas  (37),  (39),  and  (40)  yields  as  limiting  cases  the 
solution  of  the  first  boundary  value  problem  in  a wedge  of  arbitrary  angle 
[9],  or  in  an  infinite  strip  [4],  or  in  a domain  bounded  by  two  tangent  circles. 

For  the  special  values  a - -1T/4,  b*~fT/4  the  representation  reduces  to  the 
usual  one  for  the  unit  circle,  since  the  integral  (39)  can  be  evaluated  in 
closed  form  when  TT/2.  Another  special  case  of  particular  interest  is 
that  of  the  semi-circle,  for  which  a «*  0,  b-TT/4.  We  obtain,  in  general, 
a wide  variety  of  explicit  examples  for  the  study  of  the  multiple- valued 
character  of  biharmonic  functions  near  a corner  of  a curve  along  which  t^e 
boundary  conditions  (15)  hold,  and  consequently  we  can  analyze  the  viscous 
flow  around  a corner.  Finally,  since  equation  (3b)  has  constant  coefficients, 
it  can  bs  solved  by  means  of  the  Fourier  transform  in  any  infinite  strip, 
and  hence  we  can  obtain  the  solution  of  the  biharmonic  problem  in  domains 
bounded  by  logarithmic  spirals . 
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